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Abstract 

First observations of the Cabibbo suppressed decays B° — > D + K~-k + it~ and 
B~ —7- D K~tt + it~ are reported using 35 pb _1 of data collected with the LHCb 
detector. Their branching fractions are measured with respect to the corresponding 
Cabibbo favored decays, from which we obtain B(B° — > D + K~n + 7r~)/B(B° — > 
D+TT--K+TT-) = (5.9 ± 1.1 ±0.5) x 10~ 2 and B(B- -> D°K-ir + ir-)/B(B- -»■ 
D°it~it + it~) = (9.4 ± 1.3 ± 0.9) x 10 -2 , where the uncertainties are statistical and 
systematic, respectively. The B~ — > D K^ir + ir~ decay is particularly interesting, 
as it can be used in a similar way to B~ — > D°K~ to measure the CKM phase 7. 
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The standard model (SM) of particle physics provides a good description of nature up 
to the TeV scale, yet many issues remain unresolved [1], including, but not limited to, 
the hierarchy problem, the preponderance of matter over antimatter in the Universe, and 
the need to explain dark matter. One of the main objectives of the LHC is to search 
for new physics beyond the Standard Model (SM) either through direct detection or 
through interference effects in b- and c-hadron decays. In the SM, the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix j2] governs the strengths of weak charged-current interactions 
and their corresponding phases. Precise measurements on the CKM matrix parameters 
may reveal deviations from the consistency that is expected in the SM, making study 
of these decays a unique laboratory in which to search for physics beyond the standard 
model. 

The most poorly constrained of the CKM parameters is the weak phase 7 = 



promising methods of measuring this phase are through the time-independent and time- 
dependent analyses of B~ -> D°K~ El [7] and B° s ->■ Df^ E]^ respectively. 
Both approaches can be extended to higher multiplicity modes, such as B° — > D°K*°, 
B~ — > D°K~7t + 7t~ [ID] and B® — > DfK ± ir + 7i^ , which could provide a comparable level 
of sensitivity. The last two decays have not previously been observed. 

In this Letter, we report first observations of the Cabibbo-suppressed (CS) B° — > 
D + K~n + n~ and B~ ->■ D°K~n + n~ decays, where D + ->■ K~n + n + and D° -)■ K~n + , 
where charge conjugation is implied throughout this Letter. These signal decays are nor- 
malized with respect to the topologically similar Cabibbo-favored (CF) B° — > D + 7t~7t + tt~ 
and B~ — > /}°7r~7r + 7r~ decays, respectively. For brevity, we use the notation X d to refer 
to the recoiling 7r~7r + 7r~ system in the CF decays and X s for the K^n+Ti^ system in the 
CS decays. 

The analysis presented here is based on 35 pb -1 of data collected with the LHCb de- 
tector in 2010. For these measurements, the most important parts of LHCb are the vertex 
detector (VELO), the charged particle tracking system, the ring imaging Cherenkov de- 
tectors (RICH) and the trigger. The VELO is instrumental in separating particles coming 
from heavy quark decays and those emerging directly from pp interactions, by providing an 
impact parameter (IP) resolution of about 16 [im + 30 |xm/pr (transverse momentum, px 
in GeV/c). The tracking system measures charged particles' momenta with a resolution 
of a p /p ~ 0.4%(0.6%) at 5 (100) GeV/c. The RICH detectors are important to identify 
kaons and suppress the large backgrounds from pions misidentified as kaons. Events are 
selected by a two-level trigger system. The first level is hardware-based, and requires 
either a large transverse energy deposition in the calorimeter system, or a high px muon 
or pair of muons detected in the muon system. The second level, the high-level trigger, 
uses simplified versions of the offline software to reconstruct decays of b- and c-hadrons 
both inclusively and exclusively. Candidates passing the trigger selections are saved and 
used for offline analysis. A more detailed description of the LHCb detector can be found 
elsewhere [H] . In this analysis the signal and normalization modes are topologically iden- 
tical, allowing loose trigger requirements to be made with small associated uncertainty. 




Its direct measurement reaches a precision of 10° — 12° [31 H]. Two 
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In particular, we exploit the fact that 6-hadrons are produced in pairs in pp collisions, and 
include events that were triggered by the decay products of either the signal 6-hadron or 
the other 6-hadron in the event. This requirement increases the efficiency of our trigger 
selection by about 80% compared to the trigger selections requiring the signal 6-hadron 
to be responsible for triggering the event, as was done in Ref. [12]. 

The selection criteria used to reconstruct the B° — > D + tt~7t + tt~ and B~ — > D°7t~7t + tt~ 
final states are described in Ref. [12]. The Cabibbo suppression results in about a factor 
of 20 lower rate. To improve the signal-to-background ratio in the CS decay modes, 
additional selection requirements are imposed, and they are applied to both the signal 
and normalization modes. The B meson candidate is required to have > 4 GeV/c, 
IP < 60 \xm with respect to its associated primary vertex (PV), where the associated PV 
is the one having the smallest impact parameter y 2 with respect to the track. We also 
require the flight distance x 2 > 144, where the x 2 is with respect to the zero flight distance 
hypothesis, and the vertex x 2 /ndf < 5, where ndf represents the number of degrees of 
freedom in the fit. The last requirement is also applied to the vertices associated with 
X<i and X s . Three additional criteria are applied only to the CS modes. First, to remove 
the peaking backgrounds from B — > DD~ , D~ — > K~tt + tt~, we veto events where the 
invariant mass, M(X S ), is within 20 MeV/c 2 of the D s mass. Information from the 
RICH is critical to reduce background from the CF decay modes. This suppression is 
accomplished by requiring the kaon in X s to have p < 100 GeV/c (above which there is 
minimal K/n separation from the RICH), and the difference in log-likelihoods between 
the kaon and pion hypotheses to satisfy A\n£(K — 7r) > 8. The latter requirement is 
determined by optimizing N$ / y/ N$ + Nb , where we assume 100 signal events (~ 1/20 
of the CF decay yields) prior to any particle identification (PID) selection requirement, 
and the combinatorial background yield, Nb, is taken from the high B- mass sideband 
(5350-5580 MeV/c 2 ). We also make a loose PID requirement of A\n£(K — 7r) < 10 on 
the pions in X s and X&- 

Selection and trigger efficiencies are determined from simulation. Events are produced 
using pythia [13] and long-lived particles are decayed using evtgen [H] . The detector re- 
sponse is simulated with geant4 [15] . The DK~tt + tt~ final states are assumed to include 
50% DK^mO)" and 20% DK^UOO)-, with smaller contributions from DK 2 (U30)-, 
DK*(1Q80)-, DK*(892)°ir- and D 1 (2A20)K~. The resonances included in the simulation 
of the Xd system are described in Ref. [12]. The relative efficiencies, including selection 
and trigger, but not PID selection, are determined to be ^b°^d+k-tt+w- I e E -*D+-K--K+ir- = 
1.05±0.04 and eB-^D°K-n+n- /cb-->d°tv-tt+-k- — 0.942±0.036, where the uncertainties are 
statistical only. The efficiencies have a small dependence on the contributing resonances 
and their daughters' masses, and we therefore do not necessarily expect the ratios to be 
equal to unity. Moreover, the additional selections on the CS modes contribute to small 
differences between the signal and normalization modes' efficiencies. 

The PID efficiencies are determined in bins of track momentum and pseudorapidity 
(rj) using the D° daughters from D*^ 1 — > 7ffD°, D° — > K~n + calibration data, where the 
particles are identified without RICH information using the charge of the soft pion, tc s . 
The kinematics of the kaon in the X s system are taken from simulation after all offline and 
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trigger selections. Applying the PID efficiencies to the simulated decays, we determine 
the efficiencies for the kaon to pass the Aln£(A — tt) > 8 requirement to be (75.9 ± 1.5)% 
for B° D + K- Ti+Tv- and (79.2 ± 1.5)% for B~ D Q R-^ + ^-. 

Invariant mass distributions for the normalization and signal modes are shown in 
Fig. [1} Signal yields are determined through unbinned maximum likelihood fits to the sum 
of signal and several background components. The signal distributions are parametrized 
as the sum of two Gaussian functions with common means, and shape parameters, cr corc 
and / C ore that represent the width and area fraction of the narrower (core) Gaussian 
portion, and r w = cr w id e /o" C ore, which is the ratio of the wider to narrower Gaussian width. 
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Figure 1: Invariant mass distributions for (a) B° — > D + tt~7t + tt~ , (b) B~ — > D 7t~7t + tt~ , 
(c) 5° — > D + K~7Y + 7Y~ and (d) B~ — ^ Z) A"~7r + 7r~ candidates from 35 pb _1 of data for 
all selected candidates. Fits as described in the text are overlaid. 



The CF modes are first fit with f com and r w constrained to the values from simulation 
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within their uncertainties, while cr core is left as a free parameter since simulation underes- 
timates the mass resolution by ~10%. For the CF decay mode fits, the background shapes 
are the same as those described in Ref. [32]. The resulting signal shape parameters from 
the CF decay fits are then fixed in subsequent fits to the CS decay modes. For er CO rc, 
the values from the CF decay fits are scaled by width correction factors (~0.95) obtained 
from MC simulations. 

For the CS decays, invariant mass shapes of specific peaking backgrounds from other 
6-hadron decays are determined from MC simulation. The largest of these backgrounds 
comes from D ( -*^n~n + n~ decays, where one of the n~ passes the AliiC(K — n) > 8 
requirement and is misidentified as a K~ . To determine the fraction of events in which 
this occurs, we use measured PID fake rates (n faking K) obtained from D*^ 1 calibration 
data (binned in (p,rj)), and apply them to each 7r~ in Dh~'k + 'k~ simulated events. A 
decay is considered a fake if either pion has p < 100 GeV/c, and a randomly generated 
number in the interval from [0, 1] is less than that pion's determined fake rate. The 
pion's mass is then replaced by the kaon's mass, and the invariant mass of the 6-hadron is 
recomputed. The resulting spectrum is then fitted using a Crystal Ball [16] lineshape and 
its parameters are fixed in fits to data. Using this method, we find the same cross-feed rate 
of (4.4 ±0.7)% for both B° ->■ D+ii-it+it- and B~ -> D°7T-7r + 7r into B° ->■ D + K~it + it~ 
and B~ — > D K~7t + tt~ , respectively, where the uncertainty includes both statistical and 
systematic sources. A similar procedure is used to obtain the D*tt~7t + 7t~ background 
yields and shapes. The background yields are obtained by multiplying the observed CF 
signal yields in data by the cross-feed rates and the fraction of background in the region 
of the mass fit (5040-5580 MeV/c 2 ). 

We also account for backgrounds from the decays B — > DD~, D~ — > K~K + tt~, 
where the K + is misidentified as a ir + . The yields of these decays are lower, but are offset 
by a larger fake rate since the PID requirement on the particles assumed to be pions is 
significantly looser (A\nC(K— n) < 10). Using the same technique as described above, the 
fake rate is found to be (24 ±2)%. The fake yield from this source is then computed from 
the product of the measured yield of B° — > D + D~ in data (161 ± 14(stat)) and the above 
fake rate. The B~ — > D°Dj yield was not directly measured, but was determined from 
known branching fractions [17] and efficiencies from simulation. Additional uncertainty 
due to these extrapolations is included in the estimated B~ — > D°D~ background yield. 

The last sources of background, which do not contribute to the signal regions, are 
from D*K~tt + 7t~ , where the soft pion or photon from the D* is lost. The shapes of these 
low mass backgrounds are taken from the fitted D*tt~7t + tt~ shapes in the D7r~7T + 7T~ mass 
fits, and the yield ratios N(D*K~7t + 7t~)/N(DK~7t + 7t~), are constrained to be equal to 
the ratios obtained from CF mode fits with a 25% uncertainty. 

The combinatorial background is assumed to have an exponential shape. A summary 
of the signal shape parameters and the specific 6-hadron backgrounds used in the CS 
signal mode fits is given in Table [TJ 

The fitted yields are 2126±69 B° D+ti-ji+ti- and 1630±57 and B~ D°it-it + it- 
events. For the CS modes, we find 90 ± 16 B° D+R-ii+tt- and 130 ± 17 B~ -»■ 
D K~n + 'K~ signal decays. The CS decay signals have significances of 7.2 and 9.0, respec- 
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tively, calculated as a/— 21n(£ / £ max ) , where £ max and C are the fit likelihoods with the 
signal yields left free and fixed to zero, respectively. In evaluating these significances, we 
remove the constraint on N(D*K~tt + 'k~)/N(DK~'k + 7t~), which would otherwise bias the 
D*K~7t + tt~ yield toward zero and inflate Cq. Varying the signal or background shapes 
or normalizations within their uncertainties has only a minor impact on the significances. 
We therefore observe for the first time the B° — > D + K~7t + tt~ and B~ — > D°K~7t + tt~ 
decay modes. 

Table 1: Summary of parameters used in the CS mass fits. Values without uncertainties 
are fixed in the CS mode fits, and values with uncertainties are included with a Gaussian 
constraint with central values and widths as indicated. 



Parameter 


D + I<-n + n- 


D°K-n + 7i- 


Mean mass (MeV/c 2 ) 
a core (MeV/c 2 ) 


5276.3 


5276.5 


15.7 


17.5 


/core 


0.88 


0.93 


& widc/^corc 


3.32 


2.82 


N(D*7nnr) 


63 ± 10 


48 ±8 


N(D*ititit) 


47 ±9 


107 ±18 


N(DD S ) 


23 ±3 


38 ±8 


N(D*Kmr)/N(DKiTiT) 


0.62 ±0.16 


1.86 ± 0.46 



The ratios of branching fractions are given by 

B(H b H c K-7i + n-) Y cs 



x e tot> 



B(H b H c ii-it+it-) Y CF 

where V CF (V cs ) are the fitted yields in the CF (CS) decay modes, and e^ t are the 
products of the relative selection and PID efficiencies discussed previously. The results 
for the branching fractions are 



B{B° D + K'n + n-) 

B(B° D+iT-ir+n-) 

B(B~ D°K-it + it-) 

B(B- D°it-it+it-) 



(5.9 ± 1.1 ±0.5) x 10" 2 , 
(9.4 ± 1.3 ±0.9) x KT 2 , 



where the first uncertainties are statistical and the second are from the systematic sources 
discussed below. 

Most systematic uncertainties cancel in the measured ratios of branching fractions; 
only those that do not are discussed below. One source of uncertainty comes from model- 
ing of the K~7t + tt~ final state. In Ref. |12j . we compared the p and px spectra of ^ from 
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Xfc and they agreed well with simulation. We have an insufficiently large data sample 
to make such a comparison in the CS signal decay modes. The departure from unity 
of the efficiency ratios obtained from simulation are due to differences in the px spectra 
between the X<i daughters in CF decays and the X s daughters in the CS decays. These 
differences depend on the contributing resonances and the daughters' masses. We take 
the full difference of the relative efficiencies from unity (4.6% for B° and 6.1% for B~) as 
a systematic uncertainty. 

The kaon PID efficiency includes uncertainties from the limited size of the data set 
used for the efficiency determination, the limited number of events in the MC sample 
over which we average, and possible systematic effects described below. The statistical 
precision is taken as the RMS width of the kaon PID efficiency distribution obtained from 
pseudo-experiments, where in each one, the kaon PID efficiencies in each (p, rj) bin are 
fluctuated about their nominal values within their uncertainties. This contributes 1.5% 
to the overall kaon PID efficiency uncertainty. We also consider the systematic error in 
using the D* data sample to determine the PID efficiency. The procedure is tested by 
comparing the kaon PID efficiency using a MC-derived efficiency matrix with the efficiency 
obtained by directly requiring A\n£(K — 7r) > 8 on the kaon from X s in the signal MC. 
The relative difference is found to be (3.6 ± 1.9)%. We take the full difference of 3.6% as 
a potential systematic error. The total kaon PID uncertainty is 3.9%. 

The fit model uncertainty includes 3% systematic uncertainty in the yields from the 
normalization modes [12]. The uncertainties in the CS signal fits are obtained by varying 
each of the signal shape parameters within the uncertainty obtained from the CF mode 
data fits. The signal shape parameter uncertainties are 2.7% for B° and 2.5% for B~ . For 
the specific 6-hadron background shapes, we obtain the uncertainty by refitting the data 
100 times, where each fit is performed with all background shapes fluctuated within their 
covariances and subsequently fixed in the fit to data (1%). The uncertainties in the yields 
from the assumed exponential shape for the combinatorial background are estimated by 
taking the difference in yields between the nominal fit and one with a linear shape for the 
combinatorial background (2%). In total, the relative yields are uncertain by 4.5% for B° 
and 4.4% for B~. 

The limited number of MC events for determining the relative efficiencies contributes 
4.1% and 3.8% to the B° and B~ branching fraction ratio uncertainties, respectively. 
Other sources of uncertainty are negligible. In total, the uncertainties on the ratio of 
branching fractions are 8.6% for B° and 9.3% for B~ . 

We have also looked at the substructures that contribute to the CS final states. Fig- 
ure[2]shows the observed distributions of (a) K~7i + 7i~ invariant mass, (b) M(D°7r + 7i~) — 
M(D°) invariant mass difference, (c) K~n + invariant mass, and (d) 7r + 7r _ invariant mass 
for B~ — > D°K~7t + 7t~ . We show events in the B mass signal region, defined to have 
an invariant mass from 5226—5326 MeV/c 2 , and events from the high-mass sideband 
(5350—5550 MeV/c 2 ), scaled by the ratio of expected background yields in the signal 
region relative to the sideband region. An excess of events is observed predominantly 
in the low K~tt + 7t~ mass region near 1300—1400 MeV/c 2 , and the number of signal 
events decreases with increasing mass. In Fig. [21(b) there appears to be an excess of 
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~10 events in the region around 550—600 MeV/c 2 , which suggests contributions from 
-Di(2420)° or D2(2A60)° meson decays. These decays can also be used for measuring 
the weak phase 7 [15] . This yield, relative to the total, is similar to what was observed 
in B~ — > .D 7r~7r + 7r~ decays |12j . Figures [2^c) and (d) show significant enhancements 
at the K*° and p° masses, consistent with decays of excited strange states, such as 
the Ki(1270)~, Ki(1400)" and i^*(1410)~. Similar distributions are observed for the 
B° — > D + K~7t + tt~ , except that no excess of events is observed near 550—600 MeV/c 2 in 
the M(D°7t + 7t~) — M(D°) invariant mass difference. 
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Figure 2: Invariant masses within the B~ D°K~ ~7t + tt system. Shown are (a) 
M(K-it + it-) , (b) M(Dn + n-) - M(D), (c) M{K--k + ), and (d) M(tt + it-) . The points 
with error bars correspond to the signal region, and the hatched histograms represent the 
scaled sideband region. 



In summary, we report first observations of the Cabibbo-suppressed decay modes 
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B° D + K~7t + tt~ and B~ — > D^K^^-k^ and measurements of their branching frac- 
tions relative to B° — > D + 7t~7t + tt~ and B~ — > Z}°7r~7r + 7r~. The B~ — > D K~7t + tt~ decay 
is particularly interesting because, with more data, it can be used to measure the weak 
phase 7, using similar techniques as for B~ — > D°K~ and B° — > D°K*°. 
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